Semiclassical dynamics and transport of the Dirac spin 
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Semiclassical theory of spin dynamics and transport is formulated using the Dirac electron model. 
This is done by constructing a wavepacket from the positive-energy electron band, and studying its 
structure and center of mass motion. The wavepacket has a minimal size equal to the Compton 
wavelength, and has self-rotation about the average spin angular momentum, which gives rise to the 
spin magnetic moment. Geometric gauge structure in the center of mass motion provides a natural 
explanation of the spin-orbit coupling and various Yafet terms. Applications to the spin-Hall and 
spin-Nernst effects are discussed. 
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I. INTRODUCTION 

The electron has been playing a central role in mod- 
ern science and technology. It has both a fundamental 
charge and spin. With the rise of spintronics, the spin 
degree of freedom comes to the fore as it is beginning 
to be employed for data processing as well as storage. |l| 
Much has been learned on how to control the spin by 
electrical, optical as well as magnetic means. Recently, 
spin transport driven by thermal gradient has also been 
demonstrated 0, Ell- 
in this paper, we present a semiclassical theory of spin 
dynamics and transport, in order to provide an intuitive 
picture and effective calculation tool for such phenomena. 
We will focus on the Dirac model, not only because it is 
fundamental to the electron, but also it arises as effective 
theory of solid-state systems such as graphene sheet Q 
and surface of topological insulators Therefore, this 
paper can serve a dual purpose: (1) to reveal the funda- 
mental nature of the electron spin, and (2) to provide a 
simple setting for understanding spin related dynamics 
and transport phenomena in solid state systems. 

The semiclassical theory is obtained by constructing 
a wavepacket in the positive energy electron band fol- 
lowing the general framework of Culcer and Niu@. We 
find that the wavepacket has a minimal size equal to the 
Compton wavelength, and has self rotation about its av- 
erage spin, much as people imagined when the spin was 
discovered 0, HI- The self-rotation also gives rise to the 
spin magnetic moment showing the fundamental orbital 
nature of the latter. The center of mass motion has a non- 
abealian geometric gauge structure, which is shown to be 
responsible for the spin-orbit coupling as well as various 
Yafet terms. This yields a spin-dependent anomalous 
velocity under an electric field, leading to the spin Hall 
effect. It also yields a spin-dependent orbital magne- 
tization that underlies the spin Nernst effect, the spin 
dependent anomalous Nernst effect. 

The paper is organized as follows. First we construct 
the wavepacket and analyze its structure and current pro- 
file. In Sec. Ill, we discuss the magnetic moment gen- 
erated by charge circulation within the wavepacket, and 



study its coupling with a weak magnetic field. In Sec. IV, 
we derive the dynamics of the center of mass, and discuss 
the relation between spin-orbit coupling and geometric 
gauge structure. Finally, we discuss the spin Hall effect 
and spin Nernst effect in Sec. V based on the non-Abclian 
Berry curvature calculated from the Dirac theory. 



II. DIRAC ELECTRON WAVEPACKET 

When the electron spin was first discovered from the 
evidence of doublets in atomic spectra, Uhlcnbcck and 
GoudsmitQ thought it as coming from the self-rotation 
of the electron charge sphere. However, the idea was 
criticized by Lorcntz[8[, who argued that the surface of 
the sphere would have to rotate with a tangential speed 
at 137 times the speed of light to produce the accurate 
spin angular momentum. Ever since, we were left with 
no choice but to accept the spin as an abstract concept. 

In 1928, Dirac formulated the Schrodingcr equation for 
a relativistic electron. Q The Dirac equation states 
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are 4x4 matri- 
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ces defined by 2 x 2 Pauli matrices <x and identity matrix. 

The eigenenergy states are 4-component plane waves, 
with a two-fold degenerate positive energy branch, 



E(q) = 



k 2 q 2 



— e(q)mc 2 



(2) 



with hq = jmv being the relativistic momentum and 
e(q) = j(v) = (1 — v 2 / c 2 )~ 1 / 2 . There is also a two-fold 
degenerate branch of negative eigenenergy —E(q). Dirac 
assumed that these states are filled to form the vacuum. 
A hole in this negative energy branch is identified as a 
positron, the antiparticle of the electron. 

The 4-component plane-wave eigenstates are called 
Dirac spinors. They can be chosen as an orthonormal 
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set. The two spinors for the positive energy branch are 
given by 



Mq)> = 



Mq)> 




(3) 



with g-t = q x ± At q = 0, they correspond to the 
two spin eigenstates with a z = ±1. 

On the other hand, the two spinors for the negative 
energy branch are given by 
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In order to have intuitive picture of spin other than ab- 
stract operator in Dirac wave equation, we study its semi- 
classical dynamics by regarding a relativistic electron as 
a wavepacket, which contains only the positive energy 
eigenstates of the Dirac equation, 

\w) = J rfqa(q,t)e lq - r [7 ?1 (q,<)| Ul (q))+7 ?2 (q,t)| W2 (q))], 

( 5 ) 

where a(q, t) = \a\e 11 ^ describes the distribution of 
the wavepacket in momentum space. The wavepacket is 
sharply peaked at the charge center q c , and is allowed to 
have an overall phase 7(q). The probability amplitudes 
i]i and 772 describe the composition of the wavepacket 
in terms of two degenerate positive energy states with 
spin up and spin down. The normalization condition of 
the wavepacket (w\w) = 1 is satisfied if J dq\a(q,t)\ 2 = 

l,|7?l| 2 + M 2 = l. 

Now we will show that using only half of the Hilbcrt 
space, the positive energy branch, to construct the 
wavepacket results in a minimum size of the wavepacket. 
This minimum size at q = is the Compton wavelength. 
To start with, we introduce a pair of projection operators, 
V = \u\){u\\ + \u 2 )(u2\ and Q = \u 3 ){u 3 \ + |u 4 )(u 4 |. One 
can see that V projects to positive energy, V\w) = \w), 
Q projects to negative energy, and V + Q = 1. 

The mean square radius A r of the wavepacket in terms 
of the projection operators V and Q is, 

A 2 = (w\r 2 \w) - (w\r\w) 2 

= (w\r(V + Q)r\w) — {w\r\w) 2 

— (w\rVr\w) — (w\r\w) 2 + (w\rQr\w) 
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FIG. 1. (Color online) Distribution of (a) probability density 
and (b) probability current density of the wavepacket with a 
Gaussian distribution a(q, t). The length scale is in units of 
the Compton wavelength A c and the color bar is from high 
density (red) to low density (blue). The profiles of (a) and 
(b) along the z-axis are plotted in (c) and (d). 
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^VrV i s t ne mean square radius of the projected position 
operator VrV, and is a positive-dehnite quantity. The 
second term is calculated as follows : 



(u;|rQr|u;) 



A, 
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where we have used the relation between the matrix el- 
ement of position operator and velocity operator. <x = 



ry a ar\ a is the spinor-averaged spin with r\ a = 



, and 



^ c = mc 1S ^ ne Compton wavelength. 

Thus, we obtain the lower bound of the mean square 
wavepacket radius as (wlrQrlw) 1 / 2 . At q c = 0, it re- 
duces to half of the Compton wavelength. We may re- 
gard this as the minimum intrinsic radius of the electron 
wavepacket. This minimum size is a consequence of using 
only half of the Hilbert space in constructing an electron 
wavepacket and it is 137 times larger than the classi- 
cal electron radius used in Lorentz's argument. Q There- 
fore, even for the tightest possible electron wavepacket, 
the electron does not have to rotate faster than the 
speed of light. To probe the wavepacket at length scales 
smaller than the Compton wavelength, the negative en- 
ergy branch has to be involved. 

In Fig. 1, we plot the probability density, probabil- 
ity current density of a wavepacket, which are defined as 
p(r) = w^(r)w(r) and j(r) = (r)caw(r). The elec- 
tron wavepacket is spin up (in the z direction) and has a 
Gaussian distribution a(q) in momentum space with zero 
mean momentum (q c = 0). A circulating current around 
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FIG. 2. Velocity distribution v(r) (in units of c) of a rotating 
wavepacket. The distance r (in units of the Compton wave- 
length A c ) measures from the center of charge. The figure 
shows the rotating wavepacket has a rigid core with a diam- 
eter equals to the Compton wavelength. 



the spin axis is clearly seen in Fig. lb, with maxima at 
r = A c . In Fig. 2, the current density of the wavepacket 
shows a rotating velocity profile, v(r) = j(r)/p(r), much 
like that of a rigid sphere (goes linearly with the radius), 
except that beyond the edge it gradually saturates to the 
speed of light. This implies a rigid core inside the self- 
rotating wavepacket. A classical analogy of this is a uni- 
formly charged, self-rotating sphere, with a diameter of 
the Compton wavelength, which is exactly the spinning 
ball picture of Uhlenbeck and Goudsmit.Q 



III. THE SPIN MAGNETIC MOMENT 

The current circulating around the spin axis of the 
wavepacket would generate a magnetic moment M = 
^ J dr(r — r c ) x j(r) where r c = (u>|r|«;) is the cen- 
ter of the wavepacket. With some algebra, one can show 
that 

— e 

M = —(w\(r — r c ) x v\w) 

= -Y^2v*M)K al 3 x Vf) a ri a (q), (8) 

a0 

expressed in terms of the matrix element of the veloc- 
ity operator v Q( g = (u a \cxx\up), and the so-called Berry 
connection TL a p = (u a \i-J^\up) . 

After putting in the velocity operator v = col in cal- 
culation, we obtain 



M 



-eh 



2me 2 (q c ) 



q c • er 



e(<7c 



(9) 



where er = iy a crr\ a is the spinor-average spin. At q c = 0, 
it reproduces the classical result, M = — -f^er = —/ib&, 
with the Bohr magneton being fj,B = 



In the following, we will show that the magnetic mo- 
ment induced by the charge circulation is characterized 
not by the canonical angular momentum but by the spin. 

The canonical angular momentum operator is defined 
as L = mr x p = mr x jV. Unlike the momentum 
p, the canonical angular momentum is not a conserved 
quantity, dL/dt =/= 0. It is the total angular momen- 
tum J = L + S that is conserved. For a self-rotating 
Dirac wavepacket, the canonical angular momentum is 
zero (when the momentum operator p acts on the wave- 
function \w), it gives Tiq, and the matrix element q Q( g = 
implies L = (w\(r — r c ) x p\w) = 0). 

In Dirac theory, spin is represented as a 4 x 4 matrix, 

£ = i ( ^ ^ / ^ 6 °^ aui ^ c avGra g c s P m by 
calculating the expectation value of the spin operator, 



i 



2e{q c 



<r + A 



a/3 
2 Qc • cr 



e(<?c 



1 



(10) 



where T, a p = (u a \^\up). It is remarkable that the aver- 
age spin calculated from the abstract spin operator has 
the same structure (inside the square bracket of Eq.©) 
as the orbital magnetic moment obtained semiclassically. 
We can therefore relate these two quantities by 



M = 



eh 



2me(q c ) 



(11) 



where the ^-factor is 2. Note that the e in the denom- 
inator can be absorbed in the relativistic mass to form 
the relativistic Bohr magneton \ib = eh/2me(q c ). With 
q c = 0, M = -ffAtsS. 

The spin therefore can be thought of as coming from 
the charge circulation of the electron wavepacket. In fact, 
the spin is related to the mechanical angular momentum 
(the mass circulation current), L TOec /j = m{w\(r — r c ) x 
v\w) = 2hll. The g-factor of 2 is then explained by the 
fact that the mechanical angular momentum calculated 
from the mass circulating current, which is proportional 
to the charge circulating current, is twice of the spin ex- 
pectation value. In a semiconductor, the g-factor can 
deviate from 2 dramatically, [l^l The origin of the anoma- 
lous g-factor can be explained as the same way coming 
from the self rotation of electron wavepacket.pl] 

In the past, there has been a number of attempts to 
find an intuitive understanding of the spin magnetic mo- 
ment within the framework of the Dirac theory. Huang 
fl2T | suggested that it can be thought of as the cur- 
rent produced by the zitterbewegung [l3[. Ohanian [l4j 
showed that the electron spin magnetic moment origi- 
nates from a circulating flow of energy of the wave field 
based on an earlier idea of Bclinfante [15l ]. These ideas 
are similar in spirits with Uhlenbeck and Goudsmit's pic- 
ture of the spin. Here, we see that the rotating charge 
model can indeed be re-established explicitly and firmly 
within the wavepacket formulation. 
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The magnetic moment obtained above exists even in 
the absence of a magnetic field. We will show that in 
the existence of an external magnetic field, the magnetic 
moment, coming from the self- rotation of the wavepacket, 
causes an energy shift in its total energy, the Zeeman 
energy. 

First we assume the external field is weak and vary- 
ing on a length scale much larger than that of the 
wavepacket. This requirement allows us to expand the 
local Hamiltonian around the position of the charge 
center r c to the first order of the gradient correction, 
H(r c ,q, c ,t) = ff (r c ,q c ,t) + (r - r c ) • (dH/dr c ). For 
a uniform magnetic field B = V x A(r c ,i), with a 
symmetric vector potential A(r, t) = x r, we have 

H (T c ,k c ,t) = ca-hk c +/3mc 2 , where k c = q c +|A(r c ,t) 
is the kinetic momentum. The energy correction due to 
the external field is given by 

5E = (w (r - r c ) • dH/dr c 

= (w\Eij(n - r^edAj/dnlw) = -M ■ B, (12) 

We therefore observe that the Zeeman energy comes from 
the energy gradient correction and is associated with the 
M defined in Eq.©. 

When both the electric field E and magnetic field B 
are present, the total energy of a wavepacket is 

E(r c , k c ) = (w\H\w) = E (k c ) - ecj)(r c ) -MB, (13) 

where Eo(k c ) is given by Eq. (|2|) with q c —> k c , and </>(r c ) 
is the scalar potential of the electric field. 



IV. THE DYNAMICS OF THE WAVEPACKET 

Surprisingly, there arc no spin-orbit coupling in the 
wavepacket energy (Eq. (|13p ). which is expected to ap- 
pear at the first order of the electric field. One can 
quantize the semiclassical Dirac electron and show that 
the spin-orbit coupling is related to the non-canonical 
wavepacket dynamics. [Ill, [ljl [13, EH The effective La- 
grangian of a wavepacket is 



+++++++++++++++++++++ 



drj 



C = ihrf— i + hk c -7l+hk c -r c --A-r c -E(r c ,k c ). (14) 
ot c 

For a Dirac electron, the Berry connection is 72. = 

x 2 

k r x <T. 



2£(£ + l) 

From the Lagrangian, one can derive the equations of 
motion for the centers of charge position and momentum, 
correct to linear order in fields, 



eh k c 

hk r = — eE x B, 



Kk c 

em 



- I E x F 

n 



(15) 

B-F^V (16) 
emc J 



where F = (J 7 ) = r^Tf]^ ? 
is called the Berry curvature. 



\2 k c CT 



F//M 



FIG. 3. Electrons are accelerated by a static electric field in 
a parallel plate capacitor. Magnetic moment M points into 
(out of) the page, giving rise to an anomalous velocity to the 
right (left). The trajectory shift is A = A c y / eVo/2mc 2 with 
a potential difference Vo = EL. 



The equation for spin precession is given by 

d- = (\ c /e)(e/H) [B + (A c /(e + 1))E x k c ] x cr, (17) 

which agrees with the Bargmann-Michel-Telegdi 
equation [2fJ. 

When only the electric field exists (B = in Eq. (flGl) ). 
we find that the wavepacket has an anomalous velocity 
in the direction of E x F, and since T oc a at low veloc- 
ity, spin-up and spin-down electrons would have opposite 
transverse velocities (see Fig.[3j). 

Notice that the r c and k c in Eq. (14) are not a canon- 
ical pair, due to the presence of the gauge potentials 1Z. 
and A. Their connections with canonical variables r and 
p are given by (valid in weak fields), fll| | 



r c = r + 7e.(7r)+£(k c )(7r), 



hk r 



77 + - Bx 72.(tt), 



(18) 



where ir = p + |A( r), and Q a = l/2(dlZ/dk a ) ■ (TZ x 
B). This is analogous to the Peierls substitution for the 
momentum. 

Now we can re-quantize the semiclassical Dirac energy 
Eq. (fT3")) and obtain the rclativistic Pauli Hamiltonian for 
all orders of velocity [l9j . 



H(r,p) = e(ir)mc 2 — e(f)(r) + ^-cr 

e 



EXTT 

(e + l)mc 



B 

(19) 

This alternative approach is intuitive when compared to 
formal procedures of block-diagonalization, such as the 
Foldy-Wouthuysen transformation. [2l| 

In Eq. (|19p . the third term is the spin-orbit coupling 
which emerges from the first-order gradient expansion of 
the scalar potential, d<fi/dr ■ 1Z. In the literature, — elZ 
has often been called an electric dipole which couples to 
the electric field to give rise to the spin-orbit energy [22|. 
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For electrons in narrow gap semiconductors, the spin- 
orbit coupling is called a Yafet termQ. This is unfor- 
tunately artificial, because its existence depends on the 
unphysical position r which depends on the choice of the 
SU(2) gauge, instead of the true position r c . The equa- 
tions of motion based on the Pauli Hamiltonian is consis- 
tent with the Dirac theory if and only if one recognizes 
this fact. 



V. SPIN HALL EFFECT AND SPIN NERNST 
EFFECT 

The presence of the Berry curvature gives the Dirac 
electron a tiny but nonzero anomalous velocity in the 
vacuum. Similar to the electron in semiconductor, such 
a Berry curvature would lead to the spin Hall effect and 
the spin Nernst effect. The discussion below relies on 
the formulation developed previously for the Hall effect 
and the Nernst effect for spinless electrons. fl8j But their 
results are strictly applicable to the present case as long 
as the electron spin is conserved. 

For spinless electrons, the Hall conductivity is given by 



e 

T 



/(k)fi,(k), 



(20) 



where /(k) is the Fermi distribution function in equi- 
librium, Q z (k) is the Abelian Berry curvature. The 
Nernst current perpendicular to the temperature gradi- 
ent is given by j x = a xy (—V y T), and the Nernst coeffi- 
cient a xy is related to the Hall conductivity a xy via 



dE 



where a xy (E) is the Hall conductivity from all of the 
states below energy E and /i is the chemical potential. 

For electrons with spins, we need to replace the Abelian 
Berry curvature f2 2 in Eq. (f2TJ)) by the non-Abelian one 
averaged over spin, (J- z ). For a Dirac electron, in the 
limit of hk « mc, we have T z = —(X 2 ./2)a z . For a 
Dirac electron gas that is not spin-polarized, the spin- 
averaged {J- z ) is zero, even though T z itself is non-zero. 
As a result, one expects neither charge Hall effect nor 
Nernst effect. 

If the electron gas is spin polarized, then (T z ) is not 
zero and one has the anomalous Hall effect (see Eq. (|20l0 . 
At the mean time, according to Eq. (f2"Tj) . there is an 
anomalous Nernst effect. Again in the small momentum 
limit, we have 



e 2 A|n 
S 2 2 



(22) 



where n is the electron density. 

At low temperature (compared to the Fermi tempera- 
ture), the Nernst coefficient and the Hall coefficient are 



related by the Mott relation (which can be derived from 
Eq. (2D), 



7T 2 k\T da xy 
3 e dfi 



(23) 



Therefore, a xy is proportional to the density of states at 
the Fermi level, dn/dfi. 

Since the Berry curvature T z oc a z at low velocity, 
spin-up and spin-down electrons would move to oppo- 
site transverse directions. Therefore, even if the elec- 
tron gas is not spin polarized, there can still be a spin 
Hall effect (see Fig. 3). This is analogous to the emer- 
gence of the spin Hall effect in bulk (non-magnetic) 
semiconductors. 

The spin Hall conductivity is given as (valid when elec- 
tron spin is conserved), 



d 3 k 



/(k)(a z ^ z ), 



(24) 



which is approximately equal to (e/2)(A 2 /2)(n/2). Sim- 
ilarly, the spin Nernst coefficient is given by 



r 2 k%T da z xy 
3 e dji 



(25) 



and al 



h A h z V 



24tt' 2 

T. For electrons in semiconductor, 



For Fermi gas at low temperature, we have a xy 



"xy — 24 

the Berry curvature has the same structure as in vac- 

2V 2 r 1 _ 
3 VW g 

1 



uum but with different coefficient, i.e. T z 

therefore the effect can be enlarged by a fac- 



(E g +AY- 



tor 



\a z 



- l 



rj], where E g is the energy gap be- 



tween conduction band and top valence band, and A is 
energy separation between the split-off band and top va- 
lence band. For example, in GaAs with E g = 1.424cV, 
A = 0.34eV, E p = 2m e V 2 /h 2 = 22.7eV, the effect is 
enlarged by 1.3 x 10 6 times. Similar to the anomalous 
Nernst effect, such a contribution is ultimately originated 
from the Berry phase correction to (now spin-dependent) 
orbital magnetization. [l8j 



VI. CONCLUSION 

We have shown that a self-rotation picture of the 
wavepacket explains the origin of the electron spin by re- 
garding the non-relativistic electron as a wavepacket at 
the bottom of the positive energy branch of the Dirac the- 
ory. The minimum size of the wavepacket equals to the 
Compton wavelength. The magnetic moment generated 
from the circulating charge current gives the Bohr mag- 
neton in non-relativistic limit, and is responsible for the 
Zeeman energy under the external fields. The g- factor of 
2 comes from the fact that the mechanical angular mo- 
mentum from the mass circulating current is twice of the 
spin expectation value. The spin-orbit coupling emerges 
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from the first-order gradient expansion of the scalar po- 
tential and is related to the Berry connection. Finally, 
the Berry curvature plays an important role in both the 
spin Hall effect and the spin Nernst effect. Although the 
predictions of our semiclassical theory can be calculated 
from the microscopic Dirac theory, it provides not only an 
intuitive conceptual view but also a quantitatively accu- 
rate theoretical framework. The method can be directly 



transplanted to Bloch electrons in crystals, making pre- 
dictions on various thermodynamics as well as transport 
phenomena, such as spin Nernst effect discussed specifi- 
cally in this paper. 

The authors would like to thank E. I. Rashba, M. 
Stone, L. Balcnt, D. Culcer, and D. Xiao for many helpful 
discussions. 
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